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Observations were made of the optical afterglow of GRB020813 (Fox, Blake 
& Price, 2002) with the KISO observatory 1.05 m Schmidt telescope and the 
Bisei astronomical observatory 1.01 m telescope. Four-band (B, V, R, and I) 
photometric data points were obtained from 2002, August 13 10:52 to 16:46 
UT, or 0.346—0.516 days after the burst. In order to investigate the early-time 
(<1 day) evolution of the afterglow, four-band light curves were produced by 
analyzing the data taken at these two astronomical observatories, as well as 
publicly released data taken by the Magellan Baade telescope (Gladders and Hall, 
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2002c). The light curves can be approximated by a broken power law, of which 
the indices are approximately 0.46 and 1.33 before and after a break at ~0.2 
days, respectively. The optical spectral index stayed approximately constant at 



O 

~0.9 over 0.17 ~ 4.07 days after the burst. Since the temporal decay index after 
the break and the spectral index measured at that time are both consistent with 
those predicted by a spherical expansion model, the early break is unlikely to be 
a jet break, but likely to represent the end of an early bump in the light curve 



as was observed in the optical afterglow of GRB021004. 
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1. Introduction 

Afterglow of a gamma-ray burst (GRB), observed in X-ray through radio frequencies, 
can be interpreted in the fireball models, wherein a shock produced by the interaction of 
relativistic ejecta with the circumburst environment expands into the surrounding medium, 
producing broad-band synchrotron emission (Meszaros & Rees 1997; Sari, Piran, & Narayan 
1998). If the GRB is collimated into a jet, the entire jet surface becomes visible to the 
observers at some time t,-. As the jet starts to expand laterally at around t,- (Rhoads 1999), 
its sweeping area increases faster than before, leading to a stronger deceleration and hence a 
faster afterglow decay. The rapid decay of some GRB afterglows, observed recently, provide 
evidence for such jet-like or collimated ejecta (Sari, Piran, & Halpern 1999). The jet model 
would relax the energy requirements on some of the more extreme GRBs by a factor of 
several hundred (Frail et al. 2001). 

The bright and long event, GRB020813, was detected on 2002 August 13, with the 
HETE-2 spacecraft (Ricker et al. 2003). The flight localization was reported in a GCN 
Position Notice at 02:48:33 UT, 4 min 14 s after the burst trigger. The subsequent ground 
analysis of the HETE-2 data produced a refined burst location, which was reported in a 
GCN Position Notice at 05:48:35 UT, 184 minutes after the burst. The location, with 90%- 
confidence error radius of 60" (due entirely to systematic errors), is centered at a 2000 = 
19 h 46 m 38 s , <5 2000 = -19°35'16" (Villasenor et al. 2002). 

The optical afterglow was found within the 60"-radius error circle at 0.078 days (112 
minutes) after the burst, at the coordinates of a 2000 = 19 h 46 m 41. s 88, <5 2000 = -19°36'05".l 
(Fox et al., 2002). Optical spectra of the afterglow taken with the Keck observatory exhibit 
numerous absorption lines, indicating a minimum redshift of 1.254 ± 0.005 (Price et al. 
2002). Early optical light curves of the afterglow are suggested to have undergone a temporal 
break at 3.5—5 hours after the burst (Bloom et al. 2002). 

The Kiso observatory of the University of Tokyo and the Bisei Astronomical Observatory 
(BAO) have established, for the first time in Japan, the capability of multi-band follow-up 
observations of GRBs (Urata et al. 2003a). The two sites serve as valuable additions to the 
world-wide optical and infrared follow-up network, because the Japan area would otherwise 
be blank for the network. 

In order to investigate the early (within 1 day of the burst) evolution of the optical 
afterglow of GRB020813 in the flux (e.g. temporal decay and existence of a jet break) and 
the spectral slope, we have analyzed the Kiso and BAO data. In addition, we analyze the 
publicly released data of the afterglow obtained by the Magellan Baade telescope (Gladders 
and Hall, 2002c). 
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2. Observations 

We carried out follow-up observations of the optical afterglow of GRB020813 with the 
1.05 m Schmidt telescope and a 2kx2k CCD Camera at the Kiso observatory, starting at 
2002 August 13 10:52 UT (0.339 days after the burst). The field of view is 51/2 x 51/2 and 
the pixel size is 1".5 square. We performed B-, V-, R-, and /-band observations using a 
system prepared for the GRBs follow-up observations (Urata et al. 2003a). We obtained the 
multi-band data as described in talbe 1, with each dataset consisting of triple frames. 

We also performed i?-band observations using the 1.01-m telescope with a Mutoh 
CV16IIE CCD camera (Kodak KAF1602E chip) at the BAO, starting at 2002 August 13 
11:44:19 UT (0.375 days after the burst). The field of view is 7'.8 x 5/2 and the pixel size is 
0".9 square. We obtained 36 frames of i?-band data, each with 60 s exposure (table 1). 

An extensive observation of the afterglow was also performed with the Magellan Baade 
6.5 m telescope, on the nights of August 13 and 14 (0.13 - 0.94 days after the burst). These 
data were obtained in the B-, V-, R-, and /-band each for 60 s, using a TEK5 camera 
(Gladders and Hall, 2002c). The four-band images, after bias subtraction and flat-fielding, 
are made publicly available by Gladders and Halls (2002c). We have hence retrieved the 
images from the ftp site introduced by them. The data consist of 4, 13, 14, and 21 frames 
in the B-, V-, R-, and /-band, respectively. 

3. Analysis and Results 
3.1. Photometry 

We processed the Kiso and BAO data by a standard method using the NOAO IRAF. 
We used appropriate calibration data for the bias-subtraction and flat-fielding corrections. 
An example of the /-band images we obtained is shown in figure 1. Thus, the afterglow is 
clearly detected in the images of all observations listed in table 1. 

The Kiso, BAO, and Baade data have been combined with median. Flux calibrations 
among the different sites were done using the APPHOT package in IRAF, referring to the 
standard stars suggested by Henden (2002). Specifically, we utilized three of them, at the 
J2000 coordinates of (19 h 46 m 41 s .1672, -19°36'00".454), (19 h 46 m 43 s .9488, -19°36'01".166), 
and (19 h 46 m 44 s .5416, -19 35'41".662). We set the aperture size to 4 time as large as the 
FWHM of objects for each data. We summarize the results of our photometry in table I. 
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3.2. Light curves 

Figure 2 shows four-band light curves of the optical afterglow based on our measure- 
ments. We thus cover a time period of 0.17 to 0.94 days after the burst in the V-, R-, /-band, 
with the Kiso, BAO and Baade datasets. The Baade datasets densely cover the early phase 
(< 0.2 days), with an additional coverage at 0.94 days. The Kiso and BAO data, in contrast, 
constrain the light curves in the intermediate (0.35 - 0.52 days) range. 

First, we tried to fit the V-, R-, and /-band light curves by a simple power law of the 
form oc t~ a , where t is the time after the burst onset and a is a constant called decay index. 
This gave a = 1.05 with reduced x 2 {x 2 l v ) of 2.44 for the V-band, a = 1.03 with x 2 l v = 
4.12 for the /?-band, and a = 0.89 with x 2 l v — 53.6 for the /-band. Thus, none of the three 
light curves are consistent with a single power-law decay. 

Next, we tried to fit the light curves with a broken power law model expressed as 

F t = /^/**)" Ql [:L-exp(-J)]/J; (1) 
J(t,U,al,a2) = (t/Q (a2 ~ al \ 

where /*, £*, a±, and a 2 are four parameters. This functional form has no physical sig- 
nificance, but it provides a good description to the GRB990510 data, with the asymptotic 
power-law indices being a\ and a 2 at early and late times, respectively (Harrison et al. 1999). 

We have successfully fitted the above function to the V-, R-, and /-band light curves. 
For the V-b&nd, we have obtained a x = 0.46 ± 0.05, a 2 = 1.35 ± 0.01, and U = 0.22 ± 0.01 
( x 2 /z/ = 0.60 with v = 11); for the /?-band, a x = 0.46 ± 0.06, a 2 = 1.33 ± 0.01, and 
f* = 0.21±0.01 (x 2 /z/ = 0.92 with v =12); for the /-band, a x = 0.30±0.02, a 2 = 1.30±0.01, 
and t* = 0.23 ±0.01 {x 2 l v — 119 w hh v =17). Thus, the decay is relatively independent of 
color, as evidenced by nearly the same values of t* ~ 0.2 days found among the three bands. 
An initial analysis of the optical light curve suggested tj ~ 3.5 — 5 hours (Bloom et al. 2002; 
Gladders & Hall 2002b). 

We fitted a single power law to the _B-band data, because we lack data before the break. 
The obtained index is a = 1.32 ± 0.03, with x 2 l v =0.26. Thus, the slope agrees with the 
values of a 2 found in the three longer wavelengths. 

The X-ray afterglow, observed by the Chandra HETG at 0.88 - 1.78 days after the burst, 
faded in brightness according to a power law, with a decay index of 1.42 ±0.05 (Vanderspek 
et al. 2002). This is close to the optical values we measured after the break. 
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3.3. Spectral flux distributions 

We have converted the BVRI magnitudes to fluxes using the effective wavelengths and 
normalizations of Fukugita et al (1995). To remove the effects of the Galactic interstellar 
extinction, we used the reddening map of Schlegel, Finkbeiner, & Davis (1998). The Galac- 
tic reddening toward the burst is significant, E(B — V)=0.101, which implies a Galactic 
extinction of A B = 0.44, A v = 0.34, A R = 0.27, and Ai = 0.20. 

In figure 3, we plot the spectral flux distributions obtained in this way using the Kiso 
data (0.43 days after the burst) and the Baade data (0.17, 0.20 and 0.94 days after the 
burst). We fitted them with a power-law function as f(y) oc u" 13 , where f(y) is the flux 
density at frequency v and f3 is the spectral index. We have obtained (3 = 0.93 ± 0.04, 
0.82 ± 0.03, 0.93 ± 0.16, and 0.91 ± 0.07 at t= 0.17, 0.20, 0.43, and 0.94 d, respectively. The 
former two values of (3 were calculated using the V-, R-, and /-band data, while the latter 
two values utilized the four-band data. Thus, the measured values of (3 are consistent with 
being constant at (3 — 0.87 ± 0.03, because fitting them with a constant yields x 2 /u = 1.82 
with v =3. 

The optical spectral index of the GRB020813 afterglow was reported as (3 = 1.06 ±0.01 
from red side Keck spectra (5600— 9400 A) taken at about t— 0.19 d (Barth et al. 2003). 
Correcting it for the Galactic extinction quoted above, we obtain (3 ~ 0.7. Levan et al. (2002) 
reported (3 = 0.8 from the HST imaging data take at t= 4.08 d. Since these measurements 
consistently imply [3=0.7^0.9, the spectral slope (3 is suggested to have remained rather 
constant over a period of t = 0.17 to 4.08 d. 

4. Discussion 

We analyzed the GRB020813 afterglow images in four bands, observed over a period 
of 0.33~0.52 days after the burst at the Kiso and Bisei observatories. We also re-analyzed 
the publicly released four-band images taken by the Baade 6.5m telescope (Gladders & Hall 
2003c). These datasets, when combined, define relatively well-sampled light curves (in B-, 
V-, R-, and /-band) covering t= 0.17 to 0.94 d. The V-, R-, and /-band light curves are 
rather similar, and can be described by a broken power-law of a\ ~ 0.4 and «2 ~ 1-3, with 
a break point at t ~ 0.2 d. Meanwhile, the spectral index of the afterglow remained at 
(3 ~0.9. As can be seen from figure 2, the present Kiso and BAO data are essential in the 
determination of the overall light curves. 

Optical light curves of several GRBs exhibit a temporal break at t ~ 1 d or later. The 
decay index of the light curve is typically ~ 1 and ~ 2, before and after the break, respec- 
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tively. In contrast, the afterglow of GRB020813 showed an unusually early break time at 0.2 
d So far, the earliest break was observed in GRB010222, at 0.73 d (Watanabe et al. 2001). 
In addition, the decay index of a± ~ 0.46 observed from the present afterglow is smaller than 
those of other afterglows, while is similar to ct\ of other bursts. These are reminiscent 
of the afterglow of GRB021004; its light curve showed an early bump characterized by a 
brightening phase lasting until ~0.07 d, followed by a dimming with a decay index of ~ 0.2 
(Urata et al. 2003b) which is close to the present a,\. 

Usually, these temporal breaks are successfully interpreted by the jet model (Rhoads 
1999; Sari, Piran, & Halpern 1999). The jet model predicts the time dependent spectral flux 
to vary as F v (t) oc z/ _/3 t~ a , wherein the temporal decay index a and the spectral index (3 are 
both determined uniquely by the electron power law spectral index p. In order to examine 
whether the jet model is applicable to the present afterglow, we utilized the relations by Sari 
et al. (1999), and calculated the values of p and a based on the observed (3. The results are 
shown in table 2, together with the employed relations. Thus, the vales of cci = 0.3 — 0.46 we 
observed are too flat to be reproduced by the jet model under any condition, while those of 
a2 = 1.30 — 1.35 agree with the model prediction assuming that we are observing a spherical 
fireball (i.e., before the break) in the frequency range below the synchrotron cooling. 

The above discussion suggests that the behavior of the present light curves after their 
apparent break at 0.2 d is consistent with what is expected before the jet break in terms of 
the standard jet model. Therefore, the observed break is probably distinct from the usual 
jet break. If so, the electron spectral index is inferred from table 4 to be p = 2.8, which falls 
at the steepest end of the distribution of p (1.4-2.8; Panaiteschu and Kumar, 2001); it is 
close to that of GRB980519 (p ~ 2.8). If the jet break exists in this burst, it should occur 
later than t=0.94 d. We suggest that the mild break observed in the present light curves 
is reminiscent of a similar break, which was observed in the decay phase of an early bump 
in the GRB021004 afterglow. The GRB021004 afterglow showed a brightening phase over 
0.05~0.07 d after the burst based on the Kiso observation, followed by a temporal break 
around 0.2 d across which the decay index changed from 0.2 to 0.7 (Urata et al. 2003b). 
In the case of GRB021004, the optical color changed over the bump, but remained constant 
across the break like in the present light curves; the behavior may be explained by a crossing 
of the typical synchrotron frequency through the optical band (Kobayashi and Zhang, 2003). 

The authors thank Gladders & Hall for their observations as well as for making their 
data publicly available. Y. Urata acknowledges support from the Japan Society for the 
Promotion of Science (JSPS) through JSPS Research Fellowships for Young Scientists. 
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Fig. 1. — An /-band image of the GRB020813 field obtained at the Kiso observatory, with a 
300-s exposure starting at 2002 August 13 11:14 UT. The afterglow is indicated by a circle 
near the image center. 
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Fig. 2. — The B-, V-, R-, and /-band light curves based on the photometry of Kiso, Bisei 
and Baade. The smooth curves represent the broken power law model (Harrison et al. 1999) 
fitted to the V-, R-, and /-band light curves. 
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Fig. 3. — The spectral flux distributions of the optical afterglow associated with GRB020813, 
determined with the Kiso data (0.43 d) and the Baade data (0.17, 0.20, and 0.94 d). The 
straight lines give the best-fit power-law function described in the text 



Table 1: Log of Kiso and Bisei follow-up observations made on 2002 August 13. 



Start Time(UT) Filter Exposure (s) mag Site 



11:58:04 


B 


300 s x 3 


19 


.96 


± 





.29 


Kiso 


13:29:06 


B 


300 s x 3 


20 


.20 


± 





,09 


Kiso 


14:57:03 


B 


300 s x 3 


20 


.30 


± 


0. 


.18 


Kiso 


13:07:08 


V 


300 s x 3 


19 


.77 


± 





.09 


Kiso 


14:35:12 


V 


300 s x 3 


19 


.77 


± 





.25 


Kiso 


10:52:26 


R 


300 s x 3 


19 


.16 


± 





.07 


Kiso 


11:44:19 


R 


60 s x 36 


19 


.35 


± 





,13 


Bisei 


12:23:20 


R 


300 s x 3 


19 


.21 


± 





.09 


Kiso 


13:51:16 


R 


300 s x 3 


19 


.55 


± 





,06 


Kiso 


11:14:23 


I 


300 s x 3 


18. 


.47 


± 





.14 


Kiso 


12:45:17 


I 


300 s x 3 


18. 


.87 


± 


0. 


.15 


Kiso 


14:13:19 


I 


300 s x 3 


18. 


.79 


± 


0. 


.10 


Kiso 



Table 2: Values of the decay index a and the electron spectral index p, calculated from the 
measured spectral index (3 = 0.91 ± 0.07. 



Frequency* 


phase 


relation^ 


a* 


P 


Vc > v opt 


sphere 


a = 3p/2,p = 2p+l 


1.37 + 0.11 


2.81 + 0.14 


V c > "opt 


jet 


a = 2/3 + 1,?? = 2/3 + 1 


2.82 + 0.14 


2.81 + 0.14 


Vc < ^opt 


sphere 


a = 3/3/2-1/2,^ = 2/3 


0.87 + 0.11 


1.81 + 0.14 


Vc < Vopt 


jet 


a = 2/3, p = 2/3 


1.82 + 0.14 


1.81 + 0.14 



*v c is the frequency at which the spectrum breaks due to synchrotron cooling, whereas u opi is the typical 
visible light frequency. 

'''The jet model relation due to Sari et al. (1999). 
•tThe decay index calculated from the spectral index [3. 



